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Abstract 

We analyze the motion of a pressure control system described by a differential 
equation with nonlocal dissipative force. This system is composed by an oscillator, 
a membrane and a constant force. We consider the dissipative memory kernel con- 
sisting of two terms. One of them is described by the Dirac delta function which 
represents a local friction, whereas for the second one we consider two types: the 
exponential and power-law functions which represent nonlocal dissipative forces. 
For these cases, one can obtain exact solutions for the displacement and velocity. 
The long-time behaviors of these quantities are also investigated. 
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1 Introduction 



Vibration is an ubiquitous phenomenon in nature and one can perceive it in 
many macroscopic systems by the senses. In the microscopic scale this phe- 
nomenon is also present, for instance, the vibration of atoms inside of each 
molecule or vibration of molecules around their equilibrium positions in a 
solid body. In engineering the vibration is present in structural and flexible 
multibody systems. Usually these systems are accompanied by dissipation pro- 
cess. General speaking the dissipation process may be local or nonlocal. The 
nonlocal friction is related to memory effect which depends on the past motion 
of a system. This kind of friction has been used to study anomalous diffusion 
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processes which have been observed in various kinds of systems such as bacte- 
rial cytoplasm motion [1], conformational fluctuations within a single protein 
molecule [2,3] and fluorescence intermittency in single enzymes [4]. Further, it 
has been applied to structural damping systems, viscoelastic materials, seismic 
and vibration isolation (see [5,6,7,8,9,10,11] and references therein). In fact, 
nonlocal friction may be used to enhance the modelling of damping systems. 

In this work we consider a system of vibration composed by an oscillator, a 
membrane and a constant force. This simple system can be used as a pressure 
control device such as shown in Fig. 1 [12]. We note that the pressure may 
be exerted by a fluid, some mechanical mechanism or a constant force such as 
the gravitational force of the earth. In this last case the motion of the system 
should be positioned, for convenience, in the same direction of the gravita- 
tional force. However, our purpose is to employ this system to investigating 
the dissipation process by the action of membrane in the system. For elastic 
membrane one can consider the system membrane-oscillator as an oscillator. 
If the system presents nonlocal dissipation due to membrane then the mass 
coupled to the oscillator will describe different motion in comparison with that 
system with local dissipation. 

In order to model the above system we consider the following differential 
equation: 

m— + / dta (t - tx) v(h) + kx = P 9(t) , (1) 
dr Jo 

where m is the effective mass, k is the spring constant, v is the velocity, ^(t) is 
the dissipative memory kernel and 9(t) is the Heaviside function. We note that 
the dissipative memory kernel ^(t) may be a very complicated function, but 
we will restrict it for some particular cases. We will consider that the friction 
force is described by a combination of local and nonlocal forces. The dissipative 
memory kernel for a local friction is represented by the Dirac 5 function, 
whereas for the nonlocal friction we consider the exponential and power-law 
functions. For this linear system, general expression for the solution can be 
obtained. We will present the explicit solutions for the relaxation function. We 
will also analyze the long-time behavior of the solutions. 



2 Exact solutions and analysis of results 

Before analyzing the details of the above system for some specific dissipative 
memory kernel, we present general expression for the solution of Eq. (1). Eq. 
(1) is linear then it can be solved by using the Laplace transform, with the 
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initial conditions x = x(0) and vq = v(0). The application of the Laplace 
transform to Eq. (1) yields 



_ x (ms + 7) + mv + p/s 

X[S) ~ mS 2 + sl+k ' [Z) 

where 7 is the Laplace transform of the damping kernel 7(t). Now we rewrite 
Eq. (2) as follows: 



x(s) = (x (ms + 7) + mvo + pj s) G(s) (3) 
where G(t) is the relaxation function and it is the Laplace inversion of 



G(s) = r . (4) 

v ' ms 2 + s 7 y + k v ' 

For simplicity, one considers xq = v = 0, then Eq. (3) becomes 



x(s) = -G(s) . (5) 
Moreover, from Eq.(5) one can obtain v(s) which is given by 



17(a) = P G(s) . (6) 
From this last equation we arrive at 

v(t) = pG(t) (7) 

and 

x(t) =p f dtxGlh) . (8) 
Jo 

We see that the velocity is proportional to the relaxation function G(t). In 
order to obtain explicit solutions for x(t) and v(t) we should give the form of 
dissipative memory kernel ^(t). We consider three cases. 

First case. We take the dissipation process as local friction. In this case the 
dissipative memory kernel is represented by the Dirac S function and the 
friction force is proportional to v(t). The function 7 is given by j(t) = 7,j5(t) 
and its Laplace transform is 7(5) = 7,5. Substituting it into Eq. (4) we obtain 
[12] 



x(t) = — \ 1 — cos (uit) H sin (uj\t) exp (—at) \ , 7^ < Amk , (9) 
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p 



v(t) = pG(t) = sin (uit) exp (—at) , 7 5 < 4mA; , 



muoi 



V 



x(t) — j- {1 — (1 + at) exp (—at)} , 7^ = Amk , 



m 



texp (—at) , 7^ = Amk , 



x(t) 



1 - 



a 



v(t) = P G(t) 
cosh (u>2t) -+ 
v(t) = sinh (u2t) exp (—at) , 7^ > Amk , 



002 



sinh (u^t) 



exp (—at) j , 7^ > 4mA; 
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(10) 

(11) 
(12) 

(13) 

(14) 



where a = 7 (5 /2m, cji = J — — a 2 and oo 2 



or 



— . The solutions present 

m. 1 



three damped cases: Underdamped 75 < 4mA;, critical damping 7^ = 4mA; and 
overdamped 7^ > 4mA;. All the three cases have the same final equilibrium 
position which is given by Xf = p/k. Their behaviors are shown in Figs. 2 and 
3. We see that the solution (11) decays faster than other cases. This means that 
the oscillator with a critical damping moves to its final equilibrium position 
faster than other cases. Further, for the underdamped and critical damping 
cases the velocity reaches a maximum value and then it decreases to zero 
continuously. 

Second case. We take the combination of a local friction represented by the 
Dirac 5 function and exponential function for the frictional memory kernel 
which is given by j(t) = 7 <5 <5(t)+7 A e _A *. We note that the exponential damping 
model has been used to model linear elastic and viscoelastic systems (see 
[11] and references therein). The Laplace transform of 7(t) is 7(5) = 7,5 + 
7a/ (s + A). Then we obtain from Eqs. (4), (7) and (8) the following results: 



and 



x 



x(t) 



^1,3+^+271 



771 {n+j) 
rj l,2+j+2n 



/'' \ to J \ III ! 



m n=0 ni 



3=0 



J- 



\ + ^)t 

m 



n+ 00 (-m 3 Y °° f_ (T^+7 A +fc) f2 y 



m 



n=0 



nl 



3=0 



X + m) *) + XtE i^ + +i+ 2n 



m 



A + ^)t 



(15) 



(16) 



where Epj (y) is the generalized Mittag-Leffler function [13,14] defined by 



y 

E " Ay) = EWnTTy 



E tl to = to = E 



13 > 0, 5 > , 
(n + k)\y n 



dy 



t Q n!r (P(n + k) + 6)- 



(17) 



(18) 
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and T (z) is the Gamma function. 



Note that relevant information of these quantities can be obtained by analyz- 
ing their asymptotic behaviors. Therefore we use the long-time limit of the 
generalized Mittag-Leffler function [14] given by 

B -oM ~ -TrW^) ■ (19) 

From Eqs. (15), (16) and (19) yield 

/ h\+ V 

(20) 



x(t) ~ f 



7 g A + 7 A ( kXt 

1 exp 



+ + k \ 7 (5 A + 7 A + A; / 



and 



v(t) ~ / A( ^ A + 7 ^ 2 exp f- - * A * - "I , (21) 

where we consider the fact that (\ + ^ s /m)t 3> 1 and (\^ s + 7 A + fc) > 
m (A + 7^/m) 2 . We see that the displacement x(t) and velocity v(t) decay ex- 
ponentially similar to those obtained from the system with an instantaneous 
friction. Moreover, we note that the constant force does not modify the expo- 
nential asymptotic form. From Eq. (20) we can obtain the final equilibrium 
position which is given by Xf = p/k. In Fig. 4 we verify the asymptotic result 
of x(t) Eq. (15) by comparing with Eq. (20), whereas in Fig. 5 we verify the 
asymptotic result of v(t) Eq. (16) by comparing with Eq. (21). The asymptotic 
results are very close to the exact ones. 

Third case. We consider the combination of instantaneous friction represented 
by the Dirac S function and long-time frictional memory kernel given by "f(t) = 
7,5<5(t) + 7 i gt _/3 , for < /3 < 1. Then, the Laplace transform of 7(£) is 7 = 
7<5 + 7/3 s/3-1 > where 7^ = 7 /3 r(l — (3) . We note that the second term of j(t) 
describes a fractional derivative on Eq. (1) which corresponds to the Caputo 
fractional derivative [13]. The solutions for the displacement and velocity can 
be obtained from Eqs. (4), (7) and (8). The results are given by 



t 2 00 f-JL t A n 00 



.2^2- 



m 



n=0 



j=0 



rp(n+j) 
"- C/ l,3+n+(l-/3)i 



m 



(22) 



and 



m l ) \ - V ^_ 

For 7 5 = (nonlocal friction) we have 
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*(«> = £± M^V* (-?* 2 - 5 ) (24) 



m " n! ^' ^ \ m 



and 



■M^S^^fr) ' (25) 

On the other hand, for 7^ = (local friction) we have 

t 2 oo f_A t 2 



and 



x(0 = ^-£ ^r^^tn (-^ ) (26) 



= £ £ L-^J-^ ( ) . (27) 



777. ^0 ^ ! \ m 



We have numerically compared these last solutions (26) and (27) with (9)- 
(14) and they give similar results. This means that the solutions (26) and (27) 
contain the underdamped, critical damping and overdamped solutions. 

Now we analyze the asymptotic behavior of the above quantities x(t) and v(t). 
By using the long-time limit of E a ^(y) (19) we obtain 



x(t) ~ I 



■ Ypsm (Trg) r (J3) 
kirti 3 



(28) 



P7^in(7r/3)r(l + /3) 

V{t) kHt^ ' (29) 

where we consider the fact that 'y s t/m ^> 1, 7^m( 1_/3 ) > 7^ ^ and kmP > 

7/375- We can see that these asymptotic results decay as a power-law and 
they are dominated by the nonlocal dissipative force; The parameter 7 5 does 
not appears in these leading terms. From Eq. (28) we can obtain the final 
equilibrium position which is given by Xf = p/k. This means that the nonlocal 
dissipative force described by a long-time memory kernel may suppress the 
presence of a local friction in the long-time limit. Moreover, the expressions 
(28) and (29) do not depend on the mass m, then the inertial term does not 
have significant influence on the long-time behavior of the system. In Fig. 6 we 
verify the asymptotic result of x(t) obtained by Eqs. (22) and (28), whereas 
in Fig. 7 we verify the asymptotic result of v(t) Eq. (23) by comparing with 
Eq. (29). The curves converge to the same behavior in the long-time limit. 
For small values of m the oscillator does not move across the final equilibrium 
position, whereas for large values of m the oscillator can move across the final 



6 



equilibrium position. In other words, the amplitude of oscillation increases 
with the increase of mass m. This result is due to the inertial term which 
has important influence for the initial movement of the system. We note that 
the solutions (22) and (23) decay slowly due to the power-law asymptotic 
behavior. These results are in contrast to those analyzed previously. 



3 Conclusion 



In this work we have investigated the motion of a particle governed by the clas- 
sical Newtonian equation (1) under the influence of the combination of local 
and nonlocal dissipative forces, linear external force given by U (x) = —kx and 
a constant load force p. This system can be used to model a pressure control 
device. In particular, we have employed the exponential and power-law func- 
tions for the dissipative memory kernel. Exact and asymptotic solutions for the 
relaxation function G(t), x(t) and v(t) have been obtained. The asymptotic 
results have permitted us to obtain the final equilibrium position Xf = p/k. In 
the case of exponential memory kernel the asymptotic results of x(t) and v(t) 
decay exponentially similar to the system described by a local friction. For 
the power-law memory kernel we have shown the asymptotic results of x(t) 
and v(t) which decay as a power-law and the leading terms are independent 
of the parameter of local friction 7 5 and mass m; This means that the long- 
time memory friction may suppress the presence of an instantaneous friction 
in the long-time limit. We note that the system designed in Fig. 1 can be 
made experimentally. Then we hope the model described by Eq. (1) may be 
used to investigate the dynamics of membranes. If necessary other kinds of 
function for dissipative memory kernel can also be used and the solutions can 
be obtained by Eq. (4). 
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Figure Captions 



Fig. 1 - The elements of a pressure control system. 

Fig. 2 - Behavior of x(t) for different values of the mass to, in arbitrary units. 
The parameters k, p and 7 5 have the following values: k = 0.25, p — 0.1 and 
7a = 0.3. The dotted line is obtained by Eq (11) with to = 0.09. 

Fig. 3 - Behavior of v(t) for different values of the mass to, in arbitrary units. 
The parameters k, p and 7,5 have the following values: k = 0.25, p — 0.1 and 
7 tf = 0.3. The dotted line is obtained by Eq (12) with to = 0.09. 

Fig. 4 - Behavior of x(t) for different values of the mass to, in arbitrary units. 
The parameters A, k, p, 7 5 and 7 A have the following values: A = 0.1, k — 0.25, 
p = 0.1, 7 5 = 0.3 and 7 A = 0.13. The dotted line is obtained by Eq (20). 

Fig. 5 - Behavior of v(t) for different values of the mass to, in arbitrary units. 
The parameters A, k, p, 7 5 and 7 A have the following values: A = 0.1, k — 0.25, 
p = 0.1, 7 5 = 0.3 and 7 A = 0.13. The dotted line is obtained by Eq (21). 

Fig. 6 - Behavior of x(t) for different values of the mass to, in arbitrary units. 
The parameters (3, k, p, 7 5 and 7^ have the following values: (3 = 0.5, k = 0.25, 
p = 0.1, 7 5 = 0.3 and 7^ = 0.13. The dotted line is obtained by Eq (28). 

Fig. 7 - Behavior of v (t) for different values of the mass to, in arbitrary units. 
The parameters (3, k, p, 7,5 and 7^ have the following values: (3 = 0.5, k = 0.25, 
p = 0.1, 7<5 = 0.3 and 7^ = 0.13. The dotted line is obtained by Eq (29). 
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